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Abstract: We present a next-to-leading order plus parton-shower event generator for the 
production of a W boson plus two bottom quarks and a jet at hadron colliders, imple¬ 
mented in the POWHEG BOX framework. Bottom-mass effects and spin correlations of the 
decay products of the W boson are fully taken into account. The code has been automat¬ 
ically generated using the two available interfaces to MadGraph4 and GoSam, the last one 
updated to a new version. We have applied the MiNLO prescription to our Wbbj calculation, 
obtaining a finite differential cross section also in the limit of vanishing jet transverse mo¬ 
mentum. Furthermore, we have compared several key distributions for Wbbj production 
with those generated with a next-to-leading order plus parton-shower event generator for 
Wbb production, and studied their factorization- and renormalization-scale dependence. 
Finally, we have compared our results with recent experimental data from the ATLAS and 
CMS Collaborations. 
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1. Introduction 

The production of a W boson in association with two b jets at hadron collider has many 
interesting experimental and theoretical facets. On the experimental side, more interest for 
this process has been recently driven by the discovery of a light scalar particle [1, 2], whose 
characteristics point to make it a suitable candidate for being the Higgs boson responsible 
for the spontaneous symmetry breaking of the Standard Model. In this respect, Wbb + X 
is an irreducible background for HW production, with the Higgs boson decaying into b 
quarks. In addition, it is also a background to single top and top-pair production, where 
the top quark decays into a Wb pair, and to many new physics searches. 

On the theoretical side, the calculation of differential cross sections at hadronic colliders 
in the presence of massive quarks is surely more challenging than with massless partons. 
In addition, the presence of logarithmic-enhanced terms of ratios of the quark mass over 
higher scales may invalidate the perturbative expansion in the strong coupling constant a s 
(see also ref. [3] for a recent review of the subject). 
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Wbb production at next-to-leading order (NLO) in QCD has been studied for a while [4, 
5, 6, 7, 8]. All these calculations were performed in the so-called 4-flavour scheme, where 
the b quark is treated as massive (except for ref. [4], where the bottom quark is mass¬ 
less) and there is no direct contribution from the b parton-distribution function in the 
incoming hadrons. Wb + X production in the 5-flavour scheme, where the b quark is 
treated as massless and there is a contribution from the b parton-distribution function, is 
discussed, for example, in refs. [9, 10, 11]. Wbb production is also available in NLO+parton- 
shower (NLO+PS) event generators such as the POWHEG BOX [12] and MC0NL0 [8]. 

In this paper we present a NLO calculation for Wbb + 1 jet production interfaced with 
the POWHEG method [13, 14] and distributed as part of the POWHEG BOX package [15]. 
Bottom-mass effects and spin correlations of the leptonic decay products of the W boson 
have been fully taken into account. In the following we will refer to this event generator 
as Wbb j. 

The Born, real, spin- and colour-correlated Born amplitudes have been generated au¬ 
tomatically using the interface of MadGraph4 [16, 17] to the POWHEG BOX [18]. The virtual 
contribution has also been computed automatically using the interface [19] to GoSam [20, 21]. 

With the straightforward use of these two interfaces we have also generated a new 
code for Wbb production at NLO, with exact spin correlations in the decay of the W boson 
into leptons. We will refer to this event generator as Wbb. In ref. [12], Wbb production was 
interfaced with the POWHEG method, and the W decay was simulated in an approximated 
way. We made several comparisons between the generator described in ref. [12] and the 
new Wbb one, studying angular and transverse-momentum distributions of the W decay 
products, and found no sizable differences. 

The choice of factorization and renormalization scale(s) for Wbb + X production is a 
debated issue in the scientific literature: in fact, it is well known that NLO corrections 
to Wbb production are quite large [6], due to the opening of gluon-initiated channels at 
NLO. A separate paper will be needed to discuss scale dependence more thoroughly. In 
this paper, we apply the MiNLO [22] procedure to our calculation, and we leave POWHEG 
and MiNLO to choose two of the scales at which the strong coupling constant is evaluated. 
This process, in fact, starts at order a 2 EM and gets contributions up to order a 2 EM at 
next-to-leading order. The advantage of using MiNLO is twofold: 

1. First of all, we are left only with the choice of the scale(s) for the primary pro¬ 
cess, i.e. the underlying flavour and kinematic configuration after the clusterization 
operated by the MiNLO procedure*. 

2. We can provide event samples that are NLO+PS accurate for Wbbj production (this 
is by construction) but that display, at the same time, a finite differential cross section 
in the limit of the accompanying jet j becoming soft or collinear, overlapping in this 
way to the Wbb results. 

The paper is organized as follows: in sec. 2 we recall all the ingredients that are 
necessary in order to build the NLO+PS code in the POWHEG BOX, and some technical detail 

*In the MiNLO framework, by primary process we mean the process before any branching has occurred, 
i.e. H production in H + jets, and Wbb production in the present case. 
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on the change of the renormalization scheme. In sec. 3, we illustrate the modifications to 
the original MiNLO procedure for the case at hand, and the scale choice(s) for the primary 
process. In sec. 4 we study some key distributions at the NLO and Les Houches event level, 
and we compare the Wbb generator with the Wbbj+MiNLO one. We present a comparison 
with CMS and ATLAS data in sec. 5 and we give our conclusions in sec. 6. 

2. Born, real and virtual contributions 

The code for the computation of the amplitudes for Wbb and Wbbj production was gen¬ 
erated using the existing interfaces of the POWHEG BOX to MadGraph4 and GoSam [20] pre¬ 
sented in refs. [18, 19]. The one-loop amplitudes are generated with the new version 2.0 
of GoSam [21], that uses QGRAF [23], FORM [24] and SPINNEY [25] for the generation of 
the Feynman diagrams. They are then computed at running time with Ninja [26, 27], 
which is a reduction program based on the Laurent expansion of the integrand [28], and 
using OneLOop [29] for the evaluation of the scalar one-loop integrals. For unstable phase- 
space points, the reduction automatically switches to Golem95 [30], that allows to com¬ 
pute the same one-loop amplitude evaluating tensor integrals. Alternatively, the tradi¬ 
tional integrand-reduction method [31], extended to D-dimensions [32], as implemented in 
SAMURAI [33], can be used. 

We point out that a numeric calculation for Wbbj was performed in ref. [34] too, with 
the W boson treated as stable, and in ref. [26], with full spin correlations in the decay of 
the W boson. In addition, there are other automated codes that can generate the virtual 
contributions (see, for example, ref. [35]). 

In order to run the code, the user has only to select the sign of the W boson, i.e. W + 
or W ~, and its leptonic decay mode, i.e. electronic, muonic or tauonic decay, in the POWHEG 
BOX input file. Other flags to control the POWHEG BOX behavior are documented in previous 
implementations and in the Docs directory. 

2.1 The decoupling and MS schemes 

When performing a fixed-order calculation with massive quarks, one can define two consis¬ 
tent renormalization schemes that describe the same physics: the usual MS scheme, where 
all flavours are treated on equal footing, and a mixed scheme [36], that we call decoupling 
scheme, in which the nif light flavours are subtracted in the MS scheme, while the heavy- 
flavour loop is subtracted at zero momentum. In this scheme, the heavy flavour decouples 
at low energies. 

The virtual contributions generated by GoSam are computed in the decoupling scheme. 
This means that, since we are dealing with bottom-quark production, we have the correct 
virtual contributions if the strong coupling constant <a s is running with 4 light flavours, and 
if the appropriate parton distribution functions (pdfs) do not include the bottom quark in 
the evolution. 

To make contact with other results expressed in terms of the MS strong coupling 
constant, running with 5 light flavours, and with pdfs with 5 flavours, we prefer to change 
our renormalization scheme and to switch to the MS one. 
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The procedure for such a switch is well known, and was discussed in ref. [37] (see 
Appendix A for a quick review of this procedure). For Wbbj production, we need to add 
to the NLO cross section, computed in the decoupling scheme, the following two terms: 

• to the qq initial-state channel 

- 2r '£ log (€N (21) 


• to the qg initial-state channel 


2 T f 


a s 

2n 



z q9 


( 2 . 2 ) 


where B qq and B qg are the squared Born amplitude for the corresponding initial states, 
/i R and /i F are the renormalization and factorization scale, respectively, and m h is the 
bottom-quark mass. 


2.2 LO and NLO comparisons 


In this section we perform a comparison between the leading order (LO) and next-to- 
leading order results for Wbbj production, using a static renormalization and factorization 
scale. We present results for W + production at the LHC at 14 TeV. Similar results can 
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Figure 1: Transverse momentum (left) and rapidity distribution (right) of the hardest jet in Wbbj 
production. NLO and LO scale-variation bands are obtained as explained in the text. In the lower 
insert, we plot the ratio of the NLO over the LO result, computed with the central scale. The 
corresponding integration statistical error is also shown. 

be obtained for W~ production. We have set the renormalization and factorization scale 
equal to (i 

Mr = Mf = M (2-3) 

and varied \i in the range fi — {m w j 2, mn w , 2m w }. Jets and other physical parameters are 
defined as in sec. 4 and a minimum cut on the hardest jet of 25 GeV has been imposed, in 
order to have finite distributions. The total cross sections for the central scale are 


cTo = 17.3 pb, 


^nlo — 28.4 pb. 


(2.4) 


















By varying the renormalization and factorization scale we have a 46% variation for the 
LO total cross and a 29% variation at NLO, showing a significative reduction of the scale 
dependence. In fig. 1 we compare two differential distributions, i.e. the transverse momen¬ 
tum and rapidity of the hardest jet, at leading and next-to-leading order. A reduction on 
the scale dependence is clear in the two panels. 

In the rest of the paper we abandon the use of a fixed scale, since we leave to MiNLO [22] 
the choice of scales. Being interested in a shower Monte Carlo simulation, the most appro¬ 
priate scale for the evaluation of the strong coupling constant associated with the emission 
of a jet is the transverse momentum of the jet itself, provided the suitable Sudakov form 
factor is attached, as done by the MiNLO procedure. We will further illustrate our use of 
the MiNLO procedure in the next section. 

3. MiNLO 

The MiNLO procedure [22] has already been applied successfully to several production pro¬ 
cesses: H/Z/Wj [38], W/Zjj [39], HVj [19], trijet [40], VLy [41] production. In H/Z/Wj 
and HVj production, a slightly modified version of the original MiNLO formalism [38] al¬ 
lowed to reach NLO accuracy for inclusive quantities with one less jet too, i.e. H/W/Z 
and HV production, respectively. With the same modified version and with the additional 
knowledge of the fixed NNLO differential cross section [42, 43], Higgs boson [44, 45] and 
Drell-Yan [46] production could be simulated at NNLO+PS accuracy. 

Unfortunately, for Wbbj production, no such modification of the MiNLO Sudakov form 
factor is known (due to the presence of coloured particles in the final state that complicates 
the structure of the resummation formulae), and we cannot demonstrate that we can gen¬ 
erate an event sample with NLO accuracy for Wbb production too. We will limit ourselves 
to show that we can generate a Wbbj event sample with finite differential cross section 
down to the transverse momentum of the hardest jet j going to zero, and to show that 
it agrees fairly well with the cross section obtained with the Wbb event generator for Wbb 
production. 

Our last remark concerns the clusterization procedure operated by MiNLO: since no 
collinear singularities are associated with the final-state b quarks, the clusterization proce¬ 
dure is not applied to the heavy quarks. In this way, if the event gets clustered, then we 
have to deal with the kinematics of a Wbb configuration, otherwise we have a Wbbj one. 

3.1 Wbb and Wbbj+MiNLO scale choice 

Wbb and Wbbj production are multi-scale processes. The investigation of the dependence 
of the differential cross section on different scale choices is fundamental to assess the uncer¬ 
tainties on the theoretical predictions and in the comparison with the experimental data. 
In sec. 4, we limit ourselves to present results at a given dynamical scale and we compute 
scale-variation bands around it. In sec. 5, we present a few results at a different scale when 
we compare our theoretical predictions with the ATLAS data. 

Our default scale choice is the following: 
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1. In the Wbb generator, we have set the renormalization and factorization scales equal 
to 

Eb 

Mr = Mf = M = — , where E B = V§ , and s = (p w + Mb + Mb) 2 , (3.1) 

where p w , p b and p b are the momenta of the W boson, of the 6 and b quark respec¬ 
tively, at the underlying-Born level, i.e. the kinematic configuration on top of which 
the POWHEG BOX attaches the hard radiated parton with the appropriate Sudakov 
form factor. 

2. In the Wbbj+MiNLO generator, we have the freedom to set the scale(s) of the primary 
process. We have fixed this scale to be the same scale Eb of eq. (3.1), where now 
p w , p h and ps are the momenta of the W, b and 6 in the primary process if there has 
been a clusterization. If the event has not been clustered by the MiNLO procedure, 
i.e. if the underlying Born Wbbj process is not clustered by MiNLO, we take as scale 
the partonic center-of-mass energy of the event. 

With these scale choices, we have good agreement between the the Wbb and Wbbj+MiNLO 
results, as will be shown in sec. 4.2^. The fact that the scale in the Wbb generator is smaller 
than the scale of the primary process in the Wbbj+MiNLO generator is a common feature 
of the MiNLO procedure. In fact, it has already been observed in vector-boson and Higgs 
boson production plus jet [22, 39, 19], and in trijet production [40]. 

4. Results 

In this section, we present our findings for the LHC at 7 TeV. We have set the b quark 
mass at the value = 4.75 GeV and we have used the MSTW2008 [47] pdf. Clearly, any 
other pdf could have been used. 

The CKM matrix has been set to 


d s 

u ( 0.97428 0.2253 
Vq KM = c 0.2252 0.97345 

t 0.00862 0.0403 

Since the experimental data for Wbb production are 
W~ production, we do the same in our plots. 

Jets are reconstructed using the anti-£; T algorithm [48] as implemented in the Fast jet 
package [49, 50], with jet radius R — 0.7. A minimum transverse-momentum cut of 1 GeV 
has been imposed to all the jets. No cuts on 6-jets have been imposed. In all our results, the 
MiNLO procedure has always been turned on and events have been showered and hadronized 
by PYTHIA (version 6.4.25) [51], with the AMBT1 tune (call to PYTUNE(340)). 

^We leave further investigation on the choice of scales to a forthcoming paper. 


0.00347 \ 
0.041 

0.999152 ) 


(4.1) 


presented as summed over W + and 
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4.1 NLO and Les Houches event comparisons 

In this section, we compare a few interesting kinematic distributions at the NLO level and 
at the Les Houches event (LHE) level, i.e. after the first hard emission generated with the 
POWHEG method. We first compare the transverse momentum and rapidity of the first 




Figure 2: Transverse-momentum distribution of the hardest jet at NLO and LHE level. 




yil yil 

Figure 3: Rapidity distribution of the hardest jet at NLO and LHE level. 

hardest jet, that are predicted by the Wbbj generator with next-to-leading order accuracy. 

In fig. 2, we plot the transverse momentum of the hardest jet at NLO and LHE 
level. The agreement is very good over a wide kinematic range. In the right panel, the 
low pr region is illustrated: here, the LHE distribution is finite and goes to zero due to 
the Sudakov form factor coming from MiNLO, applied to the primary process, and to the 
suppression factor associated to the produced radiation, i.e. the second jet, coming from 
the POWHEG Sudakov form factor. In the NLO result, the latter Sudakov form factor is 
absent, and the result increases. 

We find very good agreement also for the rapidity of the hardest jet, y jl , shown in 
fig. 3. In the right panel, a minimum p T cut of 15 GeV on jets has been imposed. 

In figs. 4 and 5 we show the transverse-momentum distribution for the hardest and 
next-to-hardest b jets. In both cases, the agreement between the NLO and LHE result is 
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Figure 4: Transverse-momentum distribution of the hardest b jet. 
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Figure 5: Transverse-momentum distribution of the next-to-hardest b jet. 
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Figure 6: Rapidity distribution of the hardest (left panel) and next-to-hardest b jet (right panel). 


at the level of a few percent, over the entire p T range. Similar conclusions hold for their 
rapidities, y hl and 7/ b2 , illustrated in fig. 6. 

The transverse-momentum distribution of the next-to-hardest jet shown in fig. 7 is 
predicted at leading-order only. The divergent behavior at small transverse momentum in 
the NLO result is clearly visible in the right panel, where, instead, the POWHEG Sudakov 
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Figure 7: Transverse-momentum distribution of the next-to-hardest jet. 


form factor damps the LHE result. 




Figure 8: Transverse-momentum (left panel) and rapidity distribution (right panel) of the charged 
lepton from IF-boson decay. 

As last example, we plot the transverse-momentum and rapidity distribution of the 
charged lepton from VF-boson decay in fig. 8, finding again very good agreement between 
the NLO and LHE level results. 

We have compared several other kinematic distributions with different cuts and we 
have always found, when expected, excellent agreement between the NLO and the LHE 
results. 

4.2 Wbb and Wbbj+MiNLO comparisons 

In this section we compare key distributions for Wbb and Wbbj production, computed 
using the Wbb and Wbbj+MiNLO generators, respectively, in order to investigate the level of 
agreement between them, when the hardest jet j becomes soft or collinear. 

The rapidity y whh of the Wbb system computed by the Wbb generator has NLO accuracy. 
In fact, this distribution receives contributions both at the Born level and from the virtual 
and real diagrams. In fig. 9 we plot such quantity, together with the same distribution as 
predicted by Wbbj+MiNLO. 
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Figure 9: Rapidity distribution of the Wbb system computed with the Wbb and Wbbj+MiNLO 
generators, at the LHE level. In the left panel, we show the scale variation of the Wbb generator 
in red, and the central value of the Wbbj generator in green. In the right panel, we plot the scale 
variation of the Wbbj+MiNLO code in green, and the central value of the Wbb code in red. In the 
lower insert, we display the ratio between the maximum and the minimum of the band over the 
respective central value. 


The bands in the plots of this section are the envelope of the distributions obtained by 
varying the renormalization and factorization scales by a factor of 2 around the reference 
scale fi of eq. (3.1), i.e. by multiplying the factorization and the renormalization scale by 
the scale factors K F and K R , respectively, where 

(K r , K f ) = (0.5,0.5), (0.5,1), (1,0.5), (1,1), (2,1), (1,2), (2, 2). (4.2) 

These variations have been computed using the POWHEG BOX reweighting procedure, that 
recomputes the weight associated with an event in a fast way. 

In the left panel of fig. 9, we show the scale variation of the Wbb generator in red, and 
the central value of the Wbbj generator in green, at the LHE level. In the right panel, 
we plot the scale variation of the Wbbj+MiNLO code in green, and the central value of the 
Wbb code in red. In the lower insert, we display the ratio between the maximum and the 
minimum of the band over the respective central value. 

The agreement between the predictions of the Wbb and Wbbj +MiNL0 generators is within 
the scale-variation bands. We remind the reader that, if we had used the POWHEG BOX result 
for Wbbj production without MiNLO, we could have not compare these distributions, since 
the rapidity of the Wbb system would have been divergent in the limit of the accompanying 
jet becoming soft or collinear with the incoming beams. 

To better illustrate the behavior of the differential cross section in the small transverse- 
momentum region, we compare the p T of the Wbb system obtained with the two generators 
in fig. 10. In Wbb production, this distribution is predicted with leading-order accuracy 
and the POWHEG Sudakov form factor attached to the radiation makes it finite in the 
small-p T region (the Wbb system recoils against the only hard jet generated by the POWHEG 
BOX). In Wbbj production, this distribution is finite due to the presence of the POWHEG 
Sudakov form factor attached to the radiation, most likely the next-to-hardest jet, and to 
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Figure 10: Same as fig. 9 but for the transverse momentum of the Wbb system. 


the MiNLO Sudakov form factor, attached to the hardest jet accompanying the Wbb system. 
Again the agreement is very good. The finite contribution to the differential cross section 
visible in the first p^ hh bin, in Wbb production, is due to events that have not radiated at 
the LHE level. This peak is diluted away when the whole shower is completed by a Monte 



.Wbb 



Figure 11: Rapidity and transverse momentum distribution of the Wbb system, generated by the 
Wbbj+MiNLO code, and showered by PYTHIA. 

Carlo program such as PYTHIA or HERWIG [52, 53], as shown in fig. 11, right panel. In this 
figure, we plot the differential cross sections as a function of y whh and p^ hh after the shower 
has been completed by PYTHIA. No hadronization has been switched on at this level, but 
we have explicitly checked that it has a negligible effect on these distributions. 

Away from the small transverse-momentum region, the differential cross section as a 
function of p^ hh is predicted at LO by the Wbb generator, and at NLO by the Wbbj one. 
This is clearly displayed in fig. 12, where the scale variation bands for Wbb production 
reaches the 70%, while they are around 40% for Wbbj production. 

5. Comparison with ATLAS and CMS data 

In this section we compare the results that we obtained with the Wbb and Wbbj generators 
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Figure 12: Same as fig. 9 but for the transverse momentum of the Wbb system. 


with the available experimental data. 



Figure 13: Total cross section for Wbb production within the cuts of eq. (5.2). In the left panel, 
the results after the shower done by PYTHIA. In the right plot, the full showered and hadronized 
results are shown. The measured value is compared with the results from the Wbb and Wbb j -f MiNLO 
generators. 

The CMS Collaboration has measured the Wbb cross section at the LHC at 7 TeV and 
the result, reported in ref. [54], is 

ci(pp -T Wbb) x BR (W -T iiv) = 0.53 ± 0.10 pb , (5.1) 

within the following experimental cuts on jets and charged leptons 

p J T > 25 GeV, | 2 / j | < 2.4, p e T > 25 GeV, |y e | < 2.1, Ai? j5e >0.5. (5.2) 

To reconstruct jets the anti-fc T algorithm with R = 0.5 was used and only events with 
exactly two jets which passed the 6-tagging requirements were taken into account. In 
fig. 13 we compare our predictions with the measured value. In the left panel, we show 
our result after the shower done by PYTHIA, and in the right panel, the same result at the 
hadronic level. In particular, with the Wbb generator we have 

a(pp Wbb ) x BR (W ->• £u e ) = 0.50 pb (5.3) 
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and with the Wbbj+MiNLO one 


a(pp -»■ Wbb ) x BR(TV -> lv t ) = 0.49 ±g-g| pb 


(5.4) 


at the hadronic level. We find very good agreement between the cross section computed 
with the Wbb generator and the Wbb j +MiNL0 one, and both are consistent with the measured 
data. 
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Figure 14: Cross sections, within the cuts described in the text, for the ATLAS measurement of 
“1 jet”, “2 jet” and “1+2 jet” Wb production. The theoretical central results and error bands have 
been corrected for the DPI effects of tab. 1. Results are at the full shower+hadron level. 

The ATLAS Collaboration reported a measurement of W+ 6-jets (W + b + X and 
W + 66 + X) cross section at 7 TeV in ref. [55]. Candidate W + 6-jets events are required 
to have exactly one high-p T electron or muon, as well as missing transverse momentum 
consistent with a neutrino from a W boson, and one or two reconstructed jets, exactly one 
of which must be 6-tagged. Events with two or more 6-tagged jets are rejected, as are events 
with three or more jets. The details of the analysis are as follows: jets are reconstructed 
using the anti-fc T algorithm, with a radius parameter R = 0.4, and are required to have a 
transverse momentum greater than 25 GeV and absolute rapidity \y j \ < 2.1. Furthermore 
the following cuts are applied 

p e T > 25 GeV, p v T > 25 GeV, |y e | < 2.5, m+ > 60 GeV, A R j:S > 0.5. (5.5) 

In fig. 14 we plot the ATLAS results for the measured cross-sections for the “1 jet”, “2 jet” 
and “1+2 jet” fiducial regions, together with our results. In the “1 jet” bin, we show the 
cross section with only one jet, that necessarily must contain at least one 6 quark, or a 6, 
or both (clustered in a single jet). In the “2 jet” bin, we plot the cross section for events 
with two jets, only one of which is 6-tagged. In the “1+2 jet” bin, there is the sum of the 
previous two cross sections. 

With the same scale choice of sec. 3.1, our predictions for the “1+2 jet” bin are 

cwbbj ( 1+2 jet) = 2 . 93+59 P b 

a Wbb (l+2 jet) = 3.12 ++ pb (5.6) 
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correction 

1 jet 

2 jet 

1+2 jet 

DPI [pb] 

i op +°- 40 
1 - uz -0.29 

0 42 + 0 - 12 

u * oz -0.09 

1 Q 4 + 0.42 
-0.30 


Table 1: Additive corrections for double-parton interactions (DPI) as estimated by the ATLAS 
Collaboration. See tab. 7 in ref. [55]. 


at the hadron level. Since neither the Wbb nor the Wbbj+MiNLO results contain the effect 
of double-parton interactions (DPI) within the same proton, we need to correct our cross 
sections for this. The computation of this contribution is beyond the aim of the present 
paper. On the other hand, the ATLAS Collaboration has estimated such effect and has pro¬ 
vided some additive values to correct for it. We collect them in tab. 1 and after correcting 
for DPI, we get 


0-wbbj+ D pi(l+2 jet) = 4.26+Q gg pb (5.7) 

cr wbb+DPi(l+2 jet) = 4.46^Q gg pb (5-8) 

where we have estimated the total errors by simply adding in quadrature the errors from 
different sources. These results, plus the ones for “1 jet” and “2 jet” are plotted in fig. 14. 
Our predictions in eqs. (5.7) and (5.8) should be compared with the measured value 

<j(l+2 jet) = 7.1 ± 1.5 pb. (5.9) 


Although the theoretical results and the experimental data are consistent between each 
other, the level of agreement is not so good: in fact the error bands overlap only marginally. 
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Figure 15: Transverse-momentum distribution of the 6-tagged jet in the “1 jet” sample. The 
theoretical results are at the shower+hadron level. No DPI corrections included. 

The ATLAS Collaboration has also measured the p T spectrum of the 6-tagged jet in 
the “1 jet” and “2 jet” samples. We plot the measured values and our theoretical results 
in figs. 15 and 16. No DPI corrections are available for these quantities, and this partially 
explains the discrepancy between theory and data, both with the Wbb and Wbbj+MiNLO 
generator. 
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Figure 16: Transverse-momentum distribution of the 6-tagged jet in the “2 jet” sample. The 
theoretical results are at the shower-f-hadron level. No DPI corrections included. 


5.1 A different scale choice 


All the results presented so far have been computed at the scales illustrated in sec. 3.1. In 
this section, we show a few results at the same scale used by the ATLAS Collaboration 
in ref. [55], where the measured cross sections are compared with the predictions of mixed 
4- and 5-flavour NLO calculations [9, 10, 11], computed using MCFM [56]. The results are 
generated using the following central renormalization and factorization dynamical scale 

9 ■ / k\2 


2 2 2 2 i / < 

= Mr = M = m ev + (p, 


2 m 

T ) I 


+ ip' 


+ 


m\ 


+ (Pt) J 


(5.10) 


The scale-variation bands computed in ref. [55] are obtained by varying [i between a quarter 
and four times the central value, while in our results we have varied the scales independently 
as discussed in sec. 4.2. 



MCFM NLO 

LHE 

PY+SWR 

MCFM+HAD 

PY+HAD 

“1 jet” 

2 ifi+°- 78 
Z * 1D -0.59 

2 ^4+°- 27 
Z - 04 -0.41 

9 QQ +0.23 
Z -^ y -0.37 

1 QQ +0-72 
-0.54 

2 40 + 0 - 22 
z - w -0.37 

“2 jet” 

1 4Q+0-42 
- L * 4t ' 3 -0.24 

1 ^9 +0-38 
1<oz -0.30 

1 p4 +°- 31 
1 * z ^-0.25 

i o*7 +0.40 
-0.23 

-i 97 +0.31 
1>z * -0.25 


Table 2: Cross sections (in pb) for the “1 jet” and “2 jet” sample, as defined in the text, computed 
at fixed next-to-leading order with MCFM (second column), at the LHE by the P0WHEG BOX (third 
column), LHE + shower effects (fourth column), MCFM corrected for hadronization effects (fifth 
column) and LHE + shower and hadronization (last column). Scale-variation bands are also shown. 


In tab. 2 we compare the MCFM results computed by the ATLAS Collaboration with our 
results at different stages. The MCFM NLO results are as obtained by simply running the 
MCFM code, i.e. we have subtracted all the DPI corrections and undone the hadronization 
corrections applied by ATLAS. In this way, we can compare the results for the MCFM NLO 
W+lb-jet and W+lb+lj production with the LHE and LHE + shower ones. We observe 
a good level of agreement among these predictions, within the scale-variation bands. In 
addition, we notice that hadronization effects (last column of the table) have a very small 
impact on the showered results, at the level of a few percent. 
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The theoretical result for “1+2 jet” production with hadronization effects included 
(i.e. the sum of the numbers in the MCFM+HAD column), plus DPI corrections, quoted 
in ref. [55], is given by 

<t( 1+2 jet) = 4.701^5 pb. (5.11) 

This value for the cross section is in agreement, within the scale-variation band, with our 
4-flavour result obtained with the Wbbj+MiNLO generator of eq. (5.7). If instead we use the 
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Figure 17: Cross sections, within the cuts described in the text, for the ATLAS measurement of 
“1 jet”, “2 jet” and “1+2 jet” Wb production, using the scale of eq. (5.10) as central scale of the 
primary process. The theoretical central results and error bands have been corrected for the DPI 
effects of tab. 1. Results are at the full shower+hadron level. 

scale in eq. (5.10) as central scale for the primary process, we get 

<7(1+2 jet) = 5.00+93 P b > (5.12) 

with all the other results for the “1 jet” and “2 jet” sample collected in fig. 17. In this figure, 
the cross sections obtained using eq. (5.10) as scale for the primary process are displayed 
in blue, with the associated statistical and DPI error bars. The bands are obtained exactly 
as in the previous section, by varying the scales around the central one. With this scale 
choice, we have a higher degree of overlapping of the variation bands with the data, even 
if the agreement is still not perfect. 

In fig. 18, we plot the p T spectrum of the 6-tagged jet in the “1 jet” and “2 jet” samples. 
The two panels of the figure are equivalent to figs. 15 and 16, right panels, but with the 
scales of eq. (5.10). Although the blue curves are 10-20% higher than the green ones in 
figs. 15 and 16, the ratio of the theoretical predictions over the data is around the 50% 
level. We recall again that DPI corrections have not been included in these figures, and 
this partially accounts for the discrepancy between theoretical results and data. 

6. Conclusions 

The production of a VP boson in association with one or more b jets is a significant back¬ 
ground for HW production, with the Higgs boson decaying into b quarks, and to single-top 
and top-pair production in the Standard Model and to many new-physics searches. 



1 jet 2 jet 1+2 jet 
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Figure 18: Transverse-momentum distribution of the 6-tagged jet in the “1 jet” and “2 jet” 
samples, using the scale of eq. (5.10) as central scale of the primary process. The theoretical results 
are at the shower+hadron level. No DPI corrections included. 

In this paper we have presented a NLO+parton-shower event generator for Wbbj pro¬ 
duction, where bottom-quark mass effects and spin correlations of the leptonic decay prod¬ 
ucts of the W boson have been fully taken into account. The code has been automatically 
generated using the two available interfaces of MadGraph4 and GoSam to the POWHEG BOX. 

We have applied the MiNLO procedure to this process and compared several relevant 
distributions with the corresponding ones generated with the NLO+parton-shower code for 
Wbb production. We have investigated in detail the kinematic region where the transverse 
momentum of the hardest jet j in Wbbj production becomes small, and we have found 
good agreement between the two codes. 

We have shown results using a dynamical scale both for the Wbb and Wbbj generators, 
and studied their factorization- and renormalization-scale dependence. We have compared 
our results with all the experimental data collected at the LHC for Wb{b) production, 
published in the literature up to now. While we found a very good agreement with the 
Wbb cross section measured by the CMS Collaboration, the agreement with the ATLAS 
data for Wb production is less satisfactory, for both our scale choices. We point out that 
the errors on the measurements are still quite large, and more precise results expected from 
the runs at higher energies will be of help in understanding the quality of the theoretical 
predictions. 

The code for Wbbj production is available in the POWHEG BOX V2 package, under the 
User-Processes-V2/Wbbj/ directory, while the code for Wbb production is available in 
the User-Processes-V2/Wbb_dec/ folder. All the information can be found at the POWHEG 
BOX web page http://powhegbox.mib.infn.it. 
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A. The change of the renormalization scheme 


Our calculation was carried out in the mixed renormalization scheme (also called decoupling 
scheme) of ref. [36], in which the light flavours nif (4 in our process) are subtracted in the MS 
scheme, while the heavy-flavour loops are subtracted at zero momentum. In this scheme 
the heavy flavour decouples at low energy. In fact, convergent heavy-flavour loops are 
suppressed by powers of the mass of the heavy flavour. The only unsuppressed contributions 
come from divergent graphs. But those are subtracted at zero external momenta, so their 
contribution is removed by renormalization for small momenta (see also ref. [57]). 


A.l The strong coupling constant 


In order to make contact between the renormalization carried out in the decoupling and 
in the MS scheme, we need to express the strong coupling constant in the mixed 

scheme at the scale /i R , with n\{ light flavours, in terms of o^^Mr) °f the MS scheme, 
with rif light flavours. 

In the mixed scheme, charge renormalization is performed at one loop with the sub¬ 
stitution 

TV 

37r 


a s = H 2e c4 nif) (Mr) 


l-4 nlf) 


(Mr) j 


b M _ f Mr 


+ Q («s) ( J 


(A.l) 


where a | is the bare coupling constant, 


6 (n) - 

°o — 


11C A — 4 T F n 
12tt 


and 


1 1 . / A \ 

- = - + log(47r) - 7 e , 


(A.2) 


(A.3) 


where dimensional regularization in d — 4 — 2e has been used. In the MS scheme the 
renormalization prescription is 

«s = Mr 4 nf) (Mr) 11 - 4 nf) (Mr) j } ° ( a i) | » (A.4) 


where ri{ = n\f + 1. 

It is straightforward to check that the couplings in the two schemes obey the well- 
known riif-flavour and nf-flavour renormalization-group equation respectively 


d4 nif) (A%) 

d log Mr 

d«s nf) (^ R ) 

d log Mr 


-^ nif) ( 4 nif ) (^ R )) 2 + 0 (« i ), 

-b^ s) (4 nf) (/i R )) + O (al) , 


(A.5) 

(A.6) 
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that can be easily derived by imposing the independence of the bare coupling under a 
renormalization-group transformation. Up to order the solutions of eqs. (A.5) and (A.6) 
are given by 


“[““Vr) = «s nif) (^) - &o nif) "s lo S ^2 > 

m z 

(A.7) 

o^Vr) = a4 nf) M - b^aj log . 

m z 

(A.8) 

Combining eqs. (A.l) and (A.4) we have 


Mm = 4%.) -1 T * ff §+o (4 ), 

(A.9) 

which is the standard MS matching condition for flavour crossing 


<ag nif )(m) = <ag nf )(m) + O (a|) . 

(A.10) 


A.2 The parton distribution functions 


Similarly to what has been done for the strong coupling constant, we need to find the 
changes to apply to the pdfs, in order to change scheme. The parton distribution functions 
fi(x , fi F ) for n f and nif massless flavours must match at the mass m of the heavy quark [58], 
i.e. when /i F = m. More specifically, in the MS subtraction scheme, they must satisfy the 
following (see ref. [37] for more details) 


3, 

II 

7T 

for i ^ h , 

(A.ll) 

i nt \x,m ) = 0, 


(A.12) 

7* { \x,m ) = 0, 


(A.13) 


where h stands for the heavy flavour. Using the Altarelli-Parisi evolution equations to¬ 
gether with the matching conditions given in eqs. (A. 11)-(A. 13), one can easily find the 
appropriate relations between the parton densities with and n f active flavours for fi of 
the order of m. 

In fact, the Altarelli-Parisi equations for the parton densities with n f = +1 flavours 

are given by 


a4 nf V) 


( n f)( 


d log /i 2 

ides of t 
and using eq. (A.8), we get 


27r 


E 


^ d£ Ant) (x 

.. J i 




(A. 14) 


Integrating both sides of this equation between m and /x F , neglecting terms of order a 2 , 





(A.15) 


Since the heavy-quark parton distribution functions vanish at fi F = m (see eqs. (A. 12) 
and (A. 13)), we can exclude them, by putting j ^ h,h in the sum, and we can write 


/} nf) (x, Hr) = fl nf) (x,m) + 


K), 


4”'> 


27r 


(m) logA 




E. 

j+h,h 


Arif) 

z J i 


/x 


m 


if’A). 


(A.16) 
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For i — h (or i = /i), eqs. (A. 12), (A. 13) and (A. 16) yield 


fh nt \ x ’V f) 




(A.17) 


which shows that the heavy-quark pdf is of order a s . Since an equation similar to eq. (A. 16) 
holds for riif flavours 

E 

j^h,h 


<A ' 18) 


it,) = fj n “\x, m) + 


4"” 1 


2n 


(m) io g 4 

m z 


we can take the difference between eqs. (A. 16) and (A. 18), and using eqs. (A. 11) and (A. 10), 
we obtain, for i ^ h,h 

f \ ni \ x ^ f) - f) = 

= 4 nif) ( m ) , ±4_ 

2ir 8 m 2 

j^h,h 



The only splitting function that depends explicitly upon the number of light flavours is the 
gluon splitting function, so that 


P^\z) - P^\z) = r F 5(1 - z) , (A.20) 

which applied to eq. (A. 19) gives 

/x F ) - /x F ) = T f log ^ /^)(x, m). (A.21) 

The final results are then 


& i] (x,n F ) 


0(a s ) 

/j nif) 0*5 Mf) + C> («s) > for 3 4 h,h,g 


4 nf) (x,/z P ) = & li \x,i2 F ) 


2 m , til 


+ O K) 


(A.22) 
(A.23) 
(A.24) 


A.3 Summarizing 

We have now all the ingredients to translate our formulae for Wbbj production computed 
in the decoupling scheme to the standard MS scheme. At the Born level, we have contri¬ 
butions coming from two-quark initial states and a quark-gluon initial state. We can then 
schematically write the contributions to the Born cross section as 


G 


B 

qq 


G 


B 

qg 


fq nif H x i, H F ) f( nif \x2, Mf) 
4 ni d®l, Mf) 4 n A^2, Mf) 


1 3 


n 

u qq ? 


^S^Vr) 

"s^Vr)] Bqg, 


(A.25) 
(A.26) 
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where B qq and B qg are the squared Born amplitude for the qq and qg initial-state channels, 
respectively, with three powers of the strong coupling constant stripped off and put in 
front. By using eqs. (A.9), (A.23) and (A.24) we can write, up to order 


cr 


qq 


a 


qg 


, A J g A 

1 - 2 T F ^ log-A 

Z7r m 2 


1-1 
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3 

l-2T F |Mog4 

Ztt m z 
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3 

2 , ill 

[ 1 + 3 T -^ l0g SJ 


B, 


B, 


99 




l + 2T F |i(ilog4-log4 
Zn \ 5 m z m z 


B t 


qg • 


qq ’ 


(A.27) 


(A.28) 


Applying the same change of scheme to the virtual and real contributions would give 
corrections of order or higher, beyond the NLO accuracy of our calculation. In summary, 
in order to change scheme, we have to: 


• add a term 

- 2r '£ los (§) s » 

to the qq channel; 


• add a term 


to the qg channel. 


2 T — 

Ztt 




(A.29) 


(A.30) 


References 

[1] ATLAS Collaboration Collaboration, G. Aad et. al ., Observation of a new particle in the 
search for the Standard Model Higgs boson with the ATLAS detector at the LHC , Phys.Lett. 
B716 (2012) 1-29, [arXiv : 1207.7214]. 

[2] CMS Collaboration Collaboration, S. Chatrchyan et. al. , Observation of a new boson at a 
mass of 125 GeV with the CMS experiment at the LHC ' Phys.Lett. B716 (2012) 30-61, 
[arXiv: 1207.7235]. 

[3] F. Maltoni, G. Ridolfi, and M. Ubiali, b-initiated processes at the LHC: a reappraisal , JHEP 
1207 (2012) 022, [arXiv : 1203.6393]. 

[4] R. K. Ellis and S. Veseli, Strong radiative corrections to Wbb production in pp collisions , 
Phys.Rev. D60 (1999) 011501, [hep-ph/9810489]. 

[5] F. Febres Cordero, L. Reina, and D. Wackeroth, NLO QCD corrections to W boson 
production with a massive b-quark jet pair at the Tevatron pp collider , Phys.Rev. D74 (2006) 
034007, [hep-ph/0606102]. 


- 21 - 



[ 6 ] F. Febres Cordero, L. Reina, and D. Wackeroth, W- and Z-boson production with a massive 
bottom-quark pair at the Large Hadron Collider , Phys.Rev. D80 (2009) 034015, 

[arXiv: 0906.1923]. 

[7] S. Badger, J. M. Campbell, and R. Ellis, QCD corrections to the hadronic production of a 
heavy quark pair and a W-boson including decay correlations , JHEP 1103 (2011) 027, 
[arXiv: 1011.6647]. 

[ 8 ] R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, R. Pittau, et. al ., W and Zj 7 * boson 
production in association with a bottom-antibottom pair , JHEP 1109 (2011) 061, 

[arXiv: 1106.6019]. 

[9] J. Campbell, R. K. Ellis, F. Maltoni, and S. Willenbrock, Production of a W boson and two 
jets with one b-quark tag , Phys. Rev. D75 (2007) 054015, [hep-ph/0611348]. 

[10] J. Campbell, R. K. Ellis, F. Febres Cordero, F. Maltoni, L. Reina, D. Wackeroth, and 

S. Willenbrock, Associated production of a W boson and one b jet , Phys. Rev. D79 (2009) 
034023, [0809.3003]. 

[11] J. Campbell, F. Caola, F. Febres Cordero, L. Reina, and D. Wackeroth, NLO QCD 
predictions for W + 1 jet and W + 2 jet production with at least one b jet at the 7 TeV LHC , 
Phys.Rev. D 86 (2012) 034021, [arXiv: 1107.3714]. 

[ 12 ] C. Oleari and L. Reina, W^^bb production in POWHEG , JHEP 1108 (2011) 061, 

[arXiv: 1105.4488]. 

[13] P. Nason, A new method for combining NLO QCD with shower Monte Carlo algorithms , 
JHEP 11 (2004) 040, [hep-ph/0409146]. 

[14] S. Frixione, P. Nason, and C. Oleari, Matching NLO QCD computations with Parton Shower 
simulations: the POWHEG method , JHEP 11 (2007) 070, [arXiv: 0709.2092]. 

[15] S. Alioli, P. Nason, C. Oleari, and E. Re, A general framework for implementing NLO 
calculations in shower Monte Carlo programs: the POWHEG BOX , JHEP 06 (2010) 043, 
[arXiv: 1002.2581]. 

[16] T. Stelzer and W. F. Long, Automatic generation of tree level helicity amplitudes , Comput. 
Phys. Commun. 81 (1994) 357-371, [hep-ph/9401258]. 

[17] J. Alwall et. al. , MadGraph/MadEvent vj: The New Web Generation , JHEP 09 (2007) 028, 
[arXiv: 0706.2334]. 

[18] J. M. Campbell, R. K. Ellis, R. Frederix, P. Nason, C. Oleari, et. al, NLO Higgs Boson 
Production Plus One and Two Jets Using the POWHEG BOX, MadGraphj and MCFM , 
JHEP 1207 (2012) 092, [arXiv: 1202 .5475]. 

[19] G. Luisoni, P. Nason, C. Oleari, and F. Tramontano, HW ± /HZ + 0 and 1 jet at NLO with 
the POWHEG BOX interfaced to GoSam and their merging within MiNLO , JHEP 1310 
(2013) 083, [arXiv: 1306.2542]. 

[20] G. Cullen, N. Greiner, G. Heinrich, G. Luisoni, P. Mastrolia, et. al ., Automated One-Loop 
Calculations with GoSam , Eur.Phys.J. C72 ( 2012 ) 1889, [arXiv: 1111.2034]. 

[ 21 ] G. Cullen, H. van Deurzen, N. Greiner, G. Heinrich, G. Luisoni, et. al ., GO SAM-2.0: a tool 
for automated one-loop calculations within the Standard Model and beyond , Eur.Phys.J. C74 
(2014), no. 8 3001, [arXiv: 1404.7096]. 


- 22 - 



[ 22 ] K. Hamilton, P. Nason, and G. Zanderighi, MINLO: Multi-Scale Improved NLO, JHEP 1210 
(2012) 155, [arXiv: 1206.3572]. 

[23] P. Nogueira, Automatic Feynman graph generation, J.Comput.Phys. 105 (1993) 279-289. 

[24] J. Kuipers, T. Ueda, J. Vermaseren, and J. Vollinga, FORM version 4-0, 
Comput.Phys.Commun. 184 (2013) 1453-1467, [arXiv: 1203.6543]. 

[25] G. Cullen, M. Koch-Janusz, and T. Reiter, Spinney: A Form Library for Felicity Spinors, 
Comput.Phys.Commun. 182 (2011) 2368-2387, [arXiv: 1008.0803]. 

[26] H. van Deurzen, G. Luisoni, P. Mastrolia, E. Mirabella, G. Ossola, et. al., Multi-leg One-loop 
Massive Amplitudes from Integrand Reduction via Laurent Expansion, JHEP 1403 (2014) 

115, [arXiv: 1312.6678]. 

[27] T. Peraro, Ninja: Automated Integrand Reduction via Laurent Expansion for One-Loop 
Amplitudes, Comput.Phys.Commun. 185 (2014) 2771-2797, [arXiv: 1403.1229]. 

[28] P. Mastrolia, E. Mirabella, and T. Peraro, Integrand reduction of one-loop scattering 
amplitudes through Laurent series expansion, JHEP 1206 (2012) 095, [arXiv: 1203.0291]. 

[29] A. van Hameren, OneLOop: For the evaluation of one-loop scalar functions, 
Comput.Phys.Commun. 182 (2011) 2427-2438, [arXiv: 1007.4716]. 

[30] G. Cullen, J. P. Guillet, G. Heinrich, T. Kleinschmidt, E. Pilon, et. al, Golem95C: A library 
for one-loop integrals with complex masses, Comput.Phys.Commun. 182 ( 2011 ) 2276-2284, 

[arXiv: 1101.5595]. 

[31] G. Ossola, C. G. Papadopoulos, and R. Pittau, Reducing full one-loop amplitudes to scalar 
integrals at the integrand level, Nucl. Phys. B763 (2007) 147-169, [hep-ph/0609007]. 

[32] R. K. Ellis, W. Giele, and Z. Kunszt, A Numerical Unitarity Formalism for Evaluating 
One-Loop Amplitudes, JHEP 0803 (2008) 003, [arXiv: 0708.2398]. 

[33] P. Mastrolia, G. Ossola, T. Reiter, and F. Tramontano, Scattering AMplitudes from 
Unitarity-based Reduction Algorithm at the Integrand-level, JHEP 1008 (2010) 080, 

[arXiv: 1006.0710]. 

[34] L. Reina and T. Schutzmeier, Towards Wbb + j at NLO with an Automatized Approach to 
One-Loop Computations, JHEP 1209 (2012) 119, [arXiv: 1110.4438]. 

[35] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, et. al., The automated computation 
of tree-level and next-to-leading order differential cross sections, and their matching to parton 
shower simulations, JHEP 1407 (2014) 079, [arXiv: 1405.0301]. 

[36] J. C. Collins, F. Wilczek, and A. Zee, Low-Energy Manifestations of Heavy Particles: 
Application to the Neutral Current, Phys.Rev. D18 (1978) 242. 

[37] M. Cacciari, M. Greco, and P. Nason, The pr spectrum in heavy flavor hadroproduction, 
JHEP 9805 (1998) 007, [hep-ph/9803400]. 

[38] K. Hamilton, P. Nason, C. Oleari, and G. Zanderighi, Merging H/W/Z + 0 and 1 jet at NLO 
with no merging scale: a path to parton shower + NNLO matching, JHEP 1305 (2013) 082, 
[arXiv: 1212.4504]. 

[39] J. M. Campbell, R. K. Ellis, P. Nason, and G. Zanderighi, W and Z bosons in association 
with two jets using the POWHEG method, JHEP 1308 (2013) 005, [arXiv: 1303.5447]. 


-23- 



[40] A. Kardos, P. Nason, and C. Oleari, Three-jet production in POWHEG , JHEP 1404 (2014) 
043, [arXiv: 1402.4001]. 

[41] L. Barze, M. Chiesa, G. Montagna, P. Nason, O. Nicrosini, et. al ., Wy production in hadronic 
collisions using the POWHEG+MiNLO method , JHEP 1412 (2014) 039, [arXiv: 1408.5766]. 

[42] M. Grazzini, NNLO predictions for the Higgs boson signal in the H -+ WW -+ Ivlv and 
H -+ ZZ -+ 4/ decay channels , JHEP 0802 (2008) 043, [arXiv:0801.3232]. 

[43] S. Catani, L. Cieri, G. Ferrera, D. de Florian, and M. Grazzini, Vector boson production at 
hadron colliders: a fully exclusive QCD calculation at NNLO , Phys.Rev.Lett. 103 (2009) 
082001, [arXiv:0903.2120]. 

[44] K. Hamilton, P. Nason, E. Re, and G. Zanderighi, NNLOPS simulation of Higgs boson 
production , JHEP 1310 (2013) 222, [arXiv: 1309.0017]. 

[45] K. Hamilton, P. Nason, and G. Zanderighi, Finite quark-mass effects in the NNLOPS 
POWHEG+MiNLO Higgs generator , arXiv: 1501.0463. 

[46] A. Karlberg, E. Re, and G. Zanderighi, NNLOPS accurate Drell-Yan production , JHEP 1409 
(2014) 134, [arXiv: 1407.2940]. 

[47] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Parton distributions for the LHC ' 
Eur. Phys. J. C63 (2009) 189-285, [arXiv: 0901.0002]. 

[48] M. Cacciari, G. P. Salam, and G. Soyez, The anti-kr jet clustering algorithm , JHEP 04 
(2008) 063, [arXiv:0802.1189]. 

[49] M. Cacciari and G. P. Salam, Dispelling the N 3 myth for the k t jet-finder , Phys.Lett. B641 
(2006) 57-61, [hep-ph/0512210]. 

[50] M. Cacciari, G. P. Salam, and G. Soyez, Fast Jet User Manual , Eur. Phys. J. C72 (2012) 1896, 

[arXiv: 1111.6097]. 

[51] T. Sjostrand, S. Mrenna, and P. Z. Skands, PYTHIA 6. j Physics and Manual , JHEP 0605 
(2006) 026, [hep-ph/0603175]. 

[52] G. Corcella et. al ., HERWIG 6: An event generator for hadron emission reactions with 
interfering gluons (including supersymmetric processes), JHEP 01 (2001) 010, 

[hep-ph/0011363]. 

[53] G. Corcella, I. Knowles, G. Marchesini, S. Moretti, K. Odagiri, et. al ., HERWIG 6.5 release 
note , hep-ph/0210213. 

[54] CMS Collaboration Collaboration, S. Chatrchyan et. al ., Measurement of the production 
cross section for a W boson and two b jets in pp collisions at yfs = 7 TeV, Phys.Lett. B735 
(2014) 204, [arXiv: 1312.6608]. 

[55] ATLAS Collaboration Collaboration, G. Aad et. al ., Measurement of the cross-section for 
W boson production in association with b-jets in pp collisions at yfs = 7 TeV with the 
ATLAS detector , JHEP 1306 (2013) 084, [arXiv: 1302.2929]. 

[56] http://mcfm.fnal.gov. 

[57] T. Appelquist and J. Carazzone, Infrared Singularities and Massive Fields , Phys.Rev. Dll 
(1975) 2856. 

[58] J. C. Collins and W.-K. Tung, Calculating Heavy Quark Distributions , Nucl.Phys. B278 
(1986) 934. 


-24- 



